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ABSTRACT 

£3 ■ 

. Aims. The formation and evolution of cosmic structures can be probed by studying the evolution of the luminosity function 
• • of AGNs, galaxies and clusters of galaxies and of the clustering of the X-ray active Universe, compared to the IR-UV active 
^ Universe. 

, Methods. To this purpose, we have surveyed with XMM-Newton the central ~ 0.6 deg region of the ELAIS-S1 field down to 
^ ■ flux limits of ~ 5.5 x 10" 16 erg cm" 2 s" 1 (0.5-2 keV, soft band, S), ~ 2 x 10" 15 erg cm -2 s" 1 (2-10 keV, hard band, H), and 
C3 [ ~ 4 x 10" 15 erg cm" 2 s" 1 (5-10 keV,ultra hard band, HH). 

Results. We detect a total of 478 sources, 395 and 205 of which detected in the S and H bands respectively. We identified 7 
clearly extended sources and estimated their redshift through X-ray spectral fits with thermal models. In four cases the redshift 
is consistent with z=0.4. We have computed the angular correlation function of the sources in the S and H bands finding 
best fit correlation angles 0q = 5.2 ± 3.8 arcsec and 8q = 12.8 ± 7.8 arcsec, respectively. A rough estimate of the present-day 
correlation length ro can be obtained inverting the Limber equation and assuming an appropriate redshift distribution dN/dz. 
The results range between 12.8 and 9.8 h" 1 Mpc in the S band and between 17.9 and 13.4 h" 1 Mpc in the H band, with 30-40% 
statistical errors, assuming either smooth redshift distributions or redshift distributions with spikes accounting for the presence 
of significant structure at z=0.4. The relative density of the S band sources is higher near the clusters and groups at z~ 0.4 
and extends toward East and toward South/ West. This suggests that the structure is complex, with a size comparable to the 
full XMM-Newton field. Conversely, the highest relative source densities of the H band sources are located in the central-west 
region of the field. 
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1. Introduction 

One of the most challenging goals of modern cosmology 
is to understand how the structures in the Universe have 
been formed and evolved during the time. Both clusters of 

* based on observations made with XMM-Newton, an ESA 
science mission. 



galaxies and Active Galactic Nuclei (AGNs) shine power- 
fully in the X-ray band and therefore X-ray surveys are the 
most efficient tool to both trace the evolution of the cosmic 
web from low to high redshift and investigate the correla- 
tion between the formation and light-up of supermassive 
black holes (SMBH) in AGNs and the formation and evo- 
lution of galaxies. X-ray surveys are less biased against 



2 



Puccetti et al.: The XMM-Newton survey of the ELAIS-S1 field 



obscured AGN than optical UV surveys and thus reach 
higher surface densities (~ 1000 deg -2 at the flux limits 
reachable by Chandra and XMM-Newton in 50-100 ksec) 
than even faint optical AGN surveys (e.g., the COMBO 17 
survey reaches ~ 200 deg~ 2 with magnitude limit R<24, 
Wolf et al. 2004), and so they can be used to trace the 
Large Scale Structure (LSS) more efficiently than opti- 
cal AGNs. While there is today no doubt that tight links 
and feedbacks do exist between SMBH activity and galaxy 
evolution and between galaxy activity (nuclear and star- 
formation) and the larger scale environment, very little is 
known about how these links and feedbacks work in detail. 

One of the most important achievements of Chandra 
and XMM-Newton surveys concerns the redshift distri- 
bution of the hard X-ray selected sources. Indeed it has 
been realized that the accretion powered luminosity den- 
sity is dominated by obscured, low luminosity Seyfert 
like AGN at redshifts 0.5-1. On average, the activity of 
Seyfert like objects rises up to zss 1 and then decreases 
(Hasinger 2003, Fiore et al. 2003, Ueda et al. 2003, Barger 
et al 2005. lLa Franca et al. 2005|l . recalling the evolution 
of star-forming galaxies. On the other hand, QSO activ- 
ity rises steeply up to z=2-3 or even further, recalling the 
evolution of massive spheroids (Franceschini et al. 1999). 
The z=0.5-2 redshift range can then be considered as the 
"golden epoch of galaxy and AGN activity" . Both activ- 
ities are likely to be due to a) the availability of large 
masses of cold dust and gas in galaxies; b) frequent and 
efficient interactions which destabilize the gas and make 
it available for both star-formation and nuclear accretion 
(Cavaliere & Vittorini 2000; Menci et al. 2004). 

Large-area, high-sensitivity X-ray surveys can provide 
a test of the evolution of cosmic structures along the past 
light-cone and, therefore, are powerful probes of the bary- 
onic content of the Universe. The tools to address this is- 
sue are the study of the evolution of the luminosity func- 
tion of AGN, galaxies and clusters of galaxies and the 
study of the clustering of the X-ray active Universe, com- 
pared to the IR-UV active Universe. This comparison will 
allow a "direct" investigation of how different kinds of 
galaxy activity trace the cosmic web. The comparison be- 
tween AGN clustering, galaxy clustering and theoretical 
models for the evolution of the LSS can allow us to under- 
stand whether AGNs trace higher density peaks or higher 
mass haloes and how active sources are correlated with the 
environment. Deep multiwavelength surveys like GOODS 
(i.e. the Great Observatories Origins Deep Survey) sample 
a volume of the Universe too small to be representative at 
z=0.5-2, the epoch where the action is. Hydrodynamic, 
A— CDM, simulations for the formation and evolution of 
the structure in the Universe show that boxes of the order 
of 50 Mpc side can be fair samples of the Universe and, 
consequently, only at these scales all mass overdensities 
relative to the mean, will be sampled at their "true" oc- 
currence rates (see, e.g., |Springel et af~2 005). These linear 
scales translate to boxes of 2 and 1.5 degrees side at z=0.5 
and z=2, respectively (assuming a "Concordance" cosmol- 
ogy, i.e. H = 70 km s" 1 Mpc" 1 , Q M =0.3, n A = 0.7 ). 



As a first step into this direction we have surveyed with 
XMM-Newton the central ~ 0.6 deg 2 region of the ELAIS- 
Sl field, corresponding to a linear scale of ~18 Mpc at 
z=0.5. The ongoing COSMOS (Scoville et al., in prepara- 
tion) multiwavelength survey will be the next step to be 
pursued in the following several years. The ELAIS-S1 field 
(Rowan Robison et al. 2004, Oliver et al. 2000) has been 
already surveyed in the mid-IR (15 /im) by ISO (Lari et 
al. 2001), in the radio with the Australia Compact Array 
(1.4 GHz, Gruppioni et al. 1999), in the optical/NIR by 
several ESO and Australian telescopes, which provided 
deep optical (Berta et al. 2006) and NIR (Dias et al. in 
preparation, Buttery et al. in preparation) photometry in 
the B, V, R, I, z, and K bands. A campaign of optical 
spectroscopic observations is ongoing with VIMOS@VLT 
(La Franca et al. 2004). This field is also the main field of 
the Spitzer legacy program SWIRE and was observed by 
Spitzer in Dec 2004 down to flux limits of 3.7/LtJy at 3.6/im 
and of 0.15 mJy at 24^m (Lonsdale et al. 2003, 2004). 

This paper presents the XMM-Newton observations of 
the ELAIS-S1 field, the number counts in different en- 
ergy bands and the 2D clustering properties of the X-ray 
sources. The paper is organized as follows: Sect. 2 presents 
the data analysis and the X-ray source catalog; Sect. 3 
presents the number counts in four energy bands; Sect. 
4 discusses briefly the source X-ray properties; Sect. 5 
presents an analysis of the source clustering and finally 
Sect. 6 discusses our main findings. 

A H = 70 km s" 1 Mpc" 1 , Q M =0.3, fl A = 0.7 cos- 
mology is adopted throughout the paper. 

2. The X-ray observations 

A mosaic of four partially overlapping deep XMM-Newton 
pointings covers a large (~0.6 deg 2 ) and contiguous area 
of the ELAIS-S1 region. The pointings are named ELAIS- 
Sl-A (R.A.=8.91912, Dec.=-43. 31344, J2000), ELAIS-S1- 
B (R.A.=8.92154, Dec.=-43.65575, J2000), ELAIS-S1-C 
(R.A.=8.42195, Dec.=-43.30488, J2000) and ELAIS-S1-D 
(R.A.=8.42375, Dec.=-43.65327, J2000). The X-ray ob- 
servations were performed on May 2003 through July 
2003 with the European Photon Imaging Camera (EPIC: 
one PN-CCD camera (0.5-10 keV, IStrurier et al. 2001 |l 
and two MOS-CCD cameras (MOS1, MOS2, 0.3-10 keV, 
ITurner et al. 2001jl ). Table 1 gives a log of the XMM- 
Newton observations. 

2.1. Data reduction 

The data have been processed using the XMM-Newton 
Science Analysis Survey (SAS) v.5.4.1. For the fields 
ELAIS-S1-A, ELAIS-S1-B and ELAIS-S1-C, we used the 
event files linearized with a standard reduction pipeline 
(Pipeline Processing System, PPS) at the Survey Science 
Center (SSC, University of Leicester, UK). The PPS data 
of the ELAIS-S1-D field are not available in the XMM- 
Newton archive and therefore we used the raw event files 
(i.e. the Observation Data Files, ODF), which have been 
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Instrument 


Exposure (ksec) Net exposure (ksec) J Energy range 


Count rate JJ (s _i ) 


ELAIS-S1-A OBSID=0101 111 


PN 


84 45 0.5-10 keV 


4.165±0.010 


MOS1 


85 51 0.3-10 keV 


1.348±0.005 


MOS2 


85 52 0.3-10 keV 


1.352±0.005 


ELAIS-S1-B OBSID=0801, 0901, 1001, 1601 111 


PN 


92 42 0.5-10 keV 


3.794±0.012 


MOS1 


100 52 0.3-10 keV 


1.377±0.008 


MOS2 


100 53 0.3-10 keV 


1.366±0.008 


ELAIS-S1-C OBSID=1201, 1301, 1401, 2101 111 


PN 


85 51 0.5-10 keV 


4.091±0.009 


MOS1 


97 60 0.3-10 keV 


1.214±0.005 


MOS2 


97 60 0.3-10 keV 


1.183±0.005 


ELAIS-S1-D OBSID=0401 111 


PN 


88 51 0.5-10 keV 


3.921±0.009 


MOS1 


89 57 0.3-10 keV 


1.315±0.005 


MOS2 


89 62 0.3-10 keV 


1.356±0.005 



7 Net on-axis exposure time after rejection of high background periods (see Sect. 2.1); 11 Total count rate from the whole chip 
after rejection of high background periods; /7I Observation IDs. 



linearized with the XMM-SAS pipelines, EPCHAIN and 
emchain for the PN and MOS cameras respectively. 

The fields ELAIS-S1-B and ELAIS-S1-C were both 
observed at four different times. The event files of the 
first three observations (OBSID=0801, 0901, 1001) of the 
ELAIS-S1-B field and the event files of all four observa- 
tions of the ELAIS-S1-C were merged together using the 
XMM-SAS task merge, because the roll angles are simi- 
lar (within 0.03°). The fourth observation (OBSID=1601) 
of the ELAIS-S1-B has a roll angle 0.6° from the other 
three and therefore it was analyzed independently. 

Events spread at most in two contiguous pixels for PN 
(i.e. pattern=0-4) and in four contiguous pixels for MOS 
(i.e, pattern=0-12) have been selected. Event files were 
cleaned from bad pixels (hot pixels, events out of the field 
of view, etc.) and the soft proton flares. The soft pro- 
ton flares are due to solar protons with energies less than 
a few hundred of keV. The flares can produce a count 
rate up to a factor 100 greater than the mean stationary 
background count rate. They are variable during an ob- 
servation and from observation to observation, and do not 
have a predictable spectral shape or spatial distribution on 
the detector l|Read and Ponman 2003JI . In order to remove 
periods of unwanted high background level, we located 
the flares by analyzing the light curves of the count rate 
at energies higher than 10 keV, where the X-ray sources 
contribution is negligible. We rejected the time intervals 
when the count rate is higher than 1.2 counts s _1 and 
0.3 counts s _1 for the PN and MOS cameras respectively. 
These thresholds maximize the S/N of faint sources. In 
conclusion we have rejected ~ 40 — 50% of the on source 
time. 

The observations of the four XMM-Newton fields were 
brought to a common reference frame by matching the 
positions of X-ray sources with bright, point like, op- 



tical counterparts with the optical position given by 
La Franca et al. (2004)| The systematic shifts between the 
X-ray and the optical positions were of ~ 2" -4". Figure 
n shows a mosaic of the four XMM-Newton pointings 
(ELAIS-S1-A, ELAIS-S1-B, ELAIS-S1-C, ELAIS-S1-D) 
made adding up the images of the PN and MOS cam- 
eras in the energy range 0.5-10 keV. Figure [3 shows the 
relative exposure map (see Table 1). 



2.2. Source detection 

Source detection was performed on coadded (in sky coor- 
dinates) PN+MOS1-I-MOS2 images accumulated in four 
energy bands: 0.5-10 keV (full band, F), 0.5-2 keV (soft 
band, S), 2-10 keV (hard band, H), 5-10 keV (ultra hard 
band, HH). These bands have been chosen to ease the 
comparison with previous work. The EPIC effective area 
strongly decreases after 4-5 keV and therefore the 2-10 keV 
count rate is dominated by the count rate at the lower end 
of the band, unless the spectrum is extremely hard and/or 
heavily absorbed. At the flux limit of our survey we expect 
that only a minority of the 2-10 keV sources are highly ob- 
scured (see Comastri et al. 2001). In fact, the number of 
sources detected in the 2-10 keV band is similar (within 
10 %) to the number of sources detected in the 2-5 keV 
band. 

The source detections and the X-ray photometry 
were perfomed by using the PWXDetect code, devel- 
oped at INAF - Osservatorio Astronomico di Palermo, 
following Pillitteri et al. (2006). The code is derived 
from the original ROSAT code for source detection by 



Damiani et al. (1997) and allows one to combine data 
from different EPIC cameras and data taken in differ- 
ent observations, in order to achieve the deepest sensi- 
tivity. The code is based on the analysis of the wavelet 



4 Puccetti et al.: The XMM-Newton survey of the ELAIS-S1 field 




Fig. 1. Mosaic of the four XMM-Newton pointings (ELAIS-S1-A, ELAIS-S1-B, ELAIS-S1-C, ELAIS-S1-D) obtained 
by adding up the images of the PN and MOS cameras in the energy range 0.5-10 keV. The image is exposure map 
corrected. A smoothing factor with a = 1.5 pixel (or 6") was applied. Two out of the three brightest sources in the 
field are extended sources, which are presented in Section 2.5; the other one is a type 1 AGN at z=0.143 (La Franca 
et al. 2004). 



transform (WT) of the count rate image. A WT of a 
two-dimensional image is a convolution of the image with 
a "generating wavelet" kernel, which depends on posi- 
tion and length scale. In the algorithm developed by 
Damiani et al. (1997)| the generating wavelet is a "mexi- 
can hat". The length scale is a free parameter; therefore 
this method is particularly well suited for the cases in 
which the point-spread function (PSF) is varying across 
the image. It also provides robust detections of extended 
sources. 

As a first step, a preliminary background map is com- 
puted by a suitable smoothing of the image and us- 
ing also exposure maps to handle spatial exposure non- 
uniformities. Then, a local median filter is applied to re- 
compute the background at several scales, to minimize the 
effect of point sources on the background determination. 
The WT algorithm is applied at several wavelet scales. 



Spatial maxima are selected if their heights are above 
the expected background at a chosen significant level. 
Information about position, shape and count rate of the 
sources are retained for the detection scale at which each 
source has the maximum significance. The background 
is then recomputed, excluding sources found at the first 
step and WT procedure is repeated in order to detect the 
weakest sources. The final source list with their positions, 
detection scales and photometry is then built. We have 
chosen a significance level corresponding to a probability 
of 2xl0 -5 that a maximum is a Poisson fluctuations of 
the background. This corresponds to about one spurious 
source per field. To calibrate this probability we performed 
150 simulations of each field in each energy range, assum- 
ing a background equal to the total counts in each field 
(see last column of Table 1). We then applied the WT al- 
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Fig. 2. 0.5-10 kcV exposure map of the mosaic of the four XMM-Newton pointings (ELAIS-S1-A, ELAIS-S1-B, ELAIS- 
Sl-C, ELAIS-S1-D), obtained by adding to the PN exposure map, the MOS exposure maps weighted on the relative 
effective area. 



gorithm to each simulation and evaluated the significance 
level for which we have one spurious source per field. 

To evaluate count rates, we have chosen PN as the ref- 
erence detector. The scaling factor (r) between PN, MOS1 
and MOS2 depends on the relative instrument efficiency 
and source spectral shape. We evalued this factor for each 
energy band, from the PN, MOS1, MOS2 count rate ra- 
tios, predicted for a power law model with an energy index 
ole = 0.8. The count rates in the summed image (CRs U m) 
are equivalent to PN count rates according to the following 
formula: 

CRsum = CRpn + CRmOSI • fMOSl/PN + CRmOS2 • f"MOS2/PN 

where the scaling factors r MO si/PN an d 7moS2/pn are sim- 
ilar. The exposure times used to calculate the count rates 
are obtained from exposure maps, and thus they take into 
account all spatial non-uniformities due to vignetting, chip 
geometry, CCD gaps. We have checked that the source 
count rates evalued by the PWXDetect, are consistent 
with the count rates in single PN and MOS cameras. 



448 sources are detected in the F band images with 
a probability lower than 2xl0~ 5 that they are Poisson 
background fluctuations, the weakest sources having at 
least 10 net counts in each EPIC cameras (see Table 2 
for more details). The signal to noise ratio of the weakest 
detections in the coadded PN+MOS1+MOS2 images is 
~ 3 . 395, 205 and 31 sources are detected in the coadded 
S, H and HH band images respectively (see Table 3 for 
more detail). Seven of the detected sources are clearly ex- 
tended. We discuss these sources in Sect. 2.5. The quality 
of the detections and of the X-ray centroids were always 
checked interactively. We further checked the quality of 
the WT algorithm detections by comparing the distribu- 
tions of the source counts, count rates, background, signal 
to noise ratio, and probability to that obtained using the 
SAS maximum likelihood method on the XMM-Newton 
survey of the COSMOS field (Cappelluti et al., in prepa- 
ration), obtaining consistent results. 
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Table 2. Sources detected in each of the XMM- 
ELAIS-S1 field 



Energy range 


PN counts a 


MOS counts " 


N c 


ELAIS-S1-A 


0.5-10 keV 


10-450 


10-350 


152 


0.5-2 keV 


15-370 


20-250 


135 


2-10 keV 


10-110 


10-110 


73 


5-10 keV 


25-30 


20-25 


13 


ELAIS-S1-B 


0.5-10 keV 


15-5700 


15-7000 


106 


0.5-2 keV 


15-5100 


20-5800 


102 


2-10 keV 


10-600 


10-1100 


62 


5-10 keV 


10-100 


10-140 


8 


ELAIS-Sl-C 


0.5-10 keV 


15-590 


20-450 


157 


0.5-2 keV 


10-1000 


10-830 


134 


2-10 keV 


20-190 


15-140 


69 


5-10 keV 


15-40 


5-20 


11 


ELAIS-S1-D 


0.5-10 keV 


20-420 


20-300 


126 


0.5-2 keV 


10-350 


15-260 


112 


2-10 keV 


15-130 


15-80 


54 


5-10 keV 


10-35 


10-35 
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Minimum-maximum net counts of the detected sources in the 
PN camera ( a ), and MOS1+MOS2 cameras ( b ). The counts 
are not corrected for vignetting and PSF. ( c ) Number of the 
detected sources in each of the four energy bands and in each 
field. Note that the four XMM-Newton pointings are partially 
overlapping, therefore the sum of the sources in this table is 
larger than the total number of the detected sources; note also 
that the brightest source in the ELAIS-S1-B field lies close to 
a gap between two PN CCDs. 



The source catalog, including position and 
fluxes, will be electronically available from 
VizieR. It is available also at the following site: 
: / /www.mporzio. astro.it/ELAIS-Sl This site gives 
access to a multiwavelenth catalog which also includes 
optical, near infrared and mid infrared magnitudes and 
fluxes and optical spectroscopy. The optical follow-up 
campaign is currently on going and therefore the multi- 
wavelength catalog cannot be part of this publication and 
will be properly presented in a future paper (Feruglio et 
al., in preparation). 

2.3. Count rates, fluxes and sky-coverage 

Source count rates in the four bands were extracted and 
corrected for the energy dependent vignetting (as cal- 
ibrated in orbit by ILumb et al. 2003|> and PSF (using 
the analytical approximation given by Ghizzardi 200 J and 
IGhizzardi 2002|) . The PN and MOS 4 keV PSF half power 
diameter (HPD) is ~20 and ~16 arcsec, respectively, at 
small (<2 arcmin) off-axis angles; it increases only by a 
few % at off-axis angles of 10 arcmin. Even at these large 
off-axis angles, the core can be safely considered axis- 
symmetric. 



Table 3. Detected sources 



X-ray bands 


IN 


S 


395 


ri 


OAK 

zOo 


TTTT 
1 1 I 1 


ol 


F 


448 


F+S+H+HH 


27 


F+H+HH 


1 


F+S+HH 


3 


F+S+H 


140 


F+H 


33 


H only 


4 


F+S ' 


199 


F only 


45 


S only 


26 



a Number of detected sources. 



For a sample of 40 sources in a broad range of bright- 
ness and off- axis angles, we compared the count rates eval- 
uated using the WT algorithm to the count rates mea- 
sured from spectra extracted from 35—40 arcsec radius 
regions. The count rates measured from spectra were cor- 
rected for the telescope vignetting and PSF. We found 
that the count rates obtained using the WT algorithm are 
~8% and ~4% lower than the count rates measured from 
the spectra in the F band and S band respectively. We 
then corrected the count rates evaluated using the WT 
algorithm by these amounts. The weighted mean of the 
best fit spectral indices is ci£ = 0.8 ± 0.2 (ler confidence 
level) . 

Count rates were converted to fluxes using the fol- 
lowing conversion factors: 1.6xl0 -12 erg cm~ 2 s _1 (S 
band), 8.3xl0~ 12 erg cm~ 2 s" 1 (H band) and 16.5xl0" 12 
erg cm~ 2 s _1 (HH band), which are appropriate for a 
power law spectrum with energy index ag — 0.8. These 
three conversion factors are not strongly sensitive to the 
spectral shape, due to the narrow bands: for = 0.4 and 
1.0 they change by up to ~ 1% in the S band and by up to 
10% in the H band. For an absorbed power law spectrum 
with Nh = 10 22 cm -2 and olb = 1 the conversion factors 
change by < 1% in the H band and by ~ 14% in the S 
band. The conversion factor for the F band depends more 
on the spectral shape because of the wider band. For a 
power law spectrum with = 0.8±q | it changes by up 
to 30%. To minimize this error, we therefore decided to fol- 
low a different approach for the F band: we used a variable 
conversion factor which is a function of the spectral shape, 
estimated from the hardness ratios (H-S)/(H+S). The his- 
togram of the number of sources as a function of the flux 
is given in Figure El We detect sources in the 0.5-2 keV 
and 2-10 keV bands down to flux limits of ~ 5.5 x 10~ 16 
erg cm~ 2 s" 1 and ~2x 10~ 15 erg cm~ 2 s _1 respectively. 

The "Sky coverage" defines the area of the sky covered 
down to a given flux limit, as a function of the flux. Due to 
the telescope vignetting and to the increase of the size of 
the PSF in the outer regions of the detector, the sensitivity 
decreases toward the outer detector regions. 
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Fig. 3. The distribution of the detected sources as a func- 
tion of the flux in the 0.5-2 keV band (solid line, empty 
histogram), 2-10 keV band (shaded histogram) and 5-10 
keV band (black histogram). 



The minimum number of counts needed to exceed the 
Poisson fluctuations of the background, with a probability 
threshold of 2xl0~ 5 in each detector position has been 
evaluated using the background maps in each band and 
the PSF shape as a function of the off-axis angle. This 
minimum number of counts has been divided by the ex- 
posure map and thus converted to minimum detectable 
fluxes (limiting flux) using the above defined count rate- 
flux conversion factors for the S, H and HH bands. The F 
band limiting count rates were converted to fluxes using a 
conversion factor of 3.25xl0 -12 erg cm~ 2 s _1 , which is 
appropriate for a power law spectrum with cue — 0.8. The 
sky coverage at a given flux is then obtained by adding up 
the contribution of all detector regions with a given flux 
limit. Figure 0] plots the resulting sky-coverage in the four 
bands. In each band the sky coverage curve is very steep 
near the faintest detectable flux, typically decreasing by 
an order of magnitude for a change in the limiting flux of 
~ 10%. 

The main sky coverage uncertainty is due to the un- 
known spectrum of the sources near the detection limit. 
Of course, this uncertainty increases with the width of the 
band and therefore is largest in the F band. To estimate, at 
least roughly, this uncertainty, we calculated the sky cov- 
erage also for power law spectra with ole = 0.4, ole = 1.0, 
and for absorbed power law spectra with Q£ = 0.4,1.0 
and N# = 10 22 cm~ 2 , in addition to the baseline case. At 
fluxes higher than a few xl0~ 14 erg cm~ 2 s _1 the per- 
centage deviations of the sky coverage with respect to the 
baseline case are negligible for all energy bands. At a flux 
of 10 -14 erg cm -2 s _1 the deviations are ~ 75% in the F 
band, ~ 30% in the H band and negligible for the other 
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Fig. 4. The sky coverage in the four 0.5-10 keV, 0.5-2 keV, 
2-10 keV and 5-10 keV bands. The count rate-flux con- 
version factors are appropriate for a power law spectrum 
with atE = 0.8. 



bands. At a flux of 5 x 10 15 erg cm 2 s 1 the deviations 
are ~ 100% in the F band, ~ 80% in the H band and 
negligible in the S and HH bands. 

2.4. Source confusion 

Source confusion can affect at some level deep XMM- 
Newton images. At 0.5-2 keV and 2-10 keV fluxes of 
1.5 x 10~ 15 erg cm" 2 s" 1 , and 5 x 10~ 15 erg cm -2 s -1 , 
corresponding to a sky coverage of ^0.04 square degrees, 
we expect ~ 750 — 600 sources deg -2 , using the num- 
ber counts from the compilation of Moretti et al. (2003). 
The probability of finding two sources with the above 
fluxes within 12" arcsec from each other (twice the size 
of the typical detection cell) is only ~ 2 — 2.5% for the 
S and H bands, respectively. These probabilities go up 
to ~ (4 — 5)% for fluxes corresponding to the faintest de- 
tected sources in the two bands. The spatial extension and 
asymmetry in the count distribution of the sources were 
checked interactively, to understand whether same of the 
X-ray sources can actually result from the blend of two 
or more sources. We find that this may be the case for 
~ 14 sources (i.e., 3% of the whole F sample) at fluxes 
<3 x 10 -14 erg cm ~ 2 s _1 and for 2 sources at fluxes 
~ 1.4 x 10~ 13 erg cm -2 s _1 . Nine of these sources have 
been deblended and their count rates has been measured 
interactively from the images. The remaining 7 sources 
have morphologies and surface brightnesses typical of ex- 
tended sources (2 of them are bright clusters of galaxies, 
see Sect. 2.5). We conclude that source confusion does not 
significantly affect our results, especially the X-ray num- 
ber counts presented in Sect. 3. 
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Fig. 5. 0.5-10 keV contours of the seven extended sources overlayed on the R band VIMOS images (XMMES1_145, 
XMMES1.224 and XMMES1_184, XMMES1.444, XMMES1.374, XMMES1_363,XMMES1_148). The crosses indicate 
the brightest galaxy associated with each cluster. 



Table 4. Extended sources, X-ray data 



Source name RA DEC. PNT MOS" F(0.5-10keV) KT z logLx 

J2000 J2000 net counts net counts 10~ 14 c.g.s. keV erg s" 1 



XMMES 1.145 


8.44298 


-43.29217 


3228±74 


2713±68 


14.6 


2.7±0.2 


0.21±0.01 


43.3 


XMMES 1.148 


8.44469 


-43.34471 


134±16 


105±17 


0.6 


i.s±8:S 


o.27±8:ot 


42.2 


XMMES 1.184 


8.52914 


-43.37431 


203±25 


155±22 


0.8 


1.3±0.2 


oa±h 


42.8 


XMMES 1.224 


8.61410 


-43.31646 


2448±58 


1865±51 


14.3 


3.2±0.3 


0.39±0.02 


43.9 


XMMES 1.363 


8.94811 


-43.38005 


303±25 


261±23 


1.6 


2.7±8:I 


1.10±0.13 


44.1 


XMMES 1.374 


8.97218 


-43.35345 


283±24 


224±22 


1.2 


l.liO.l 


0.37±0.06 


42.7 


XMMES 1.444 


9.11833 


-43.47706 


201±22 


176±20 


2.5 


3.i±?:f 


0.7±U 


43.7 



"Net PN counts in the 0.5-10 keV energy range. b Net MOS1+MOS2 counts in the 0.3-10 keV energy range. The quoted errors 
corresponds to 68% confidence level. 



2.5. Extended sources 

We identified seven extended X-ray sources by comparing 
their background subtracted surface brightness, computed 
in concentric annuli, to the background subtracted surface 



brightness of point-like sources found at similar off-axis 
angles, which provide a reasonably good representation 
of the local PSF. The surface brightness radial profiles of 
these sources, normalized to the peak surface brightness 
are shown in Figure For all of them, the difference with 
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Table 5. Extended sources: optical and near infrared data 1 of the brightest galaxy associated with each 
cluster. 



Source name 


R.A. 


DEC. 


B 


V 


R 


K 


J 


Zphot 




J2000 


J2000 


mag 


mag 


mag 


mag 


mag 




XMMES1.145 


8.4433322 


-43.291992 


19.24 


17.76 


16.99 


13.86 


15.15 


0.15-0.30 


XMMES 1.148 


8.4452075 


-43.344633 


20.07 


18.54 


17.74 


14.57 


15.90 


0.15-0.75 


XMMES1.184 


8.5332721 


-43.377095 


20.79 


19.31 


18.25 


14.59 


16.13 


0.25-0.80 


XMMES1.224 


8.6080313 


-43.314348 


20.79 


19.19 


17.93 


14.36 


14.55 


0.25-0.85 


XMMES1.363 


8.948075 


-43.379401 


23.31 


21.34 


20.21 


16.13 


17.48 


0.65-0.95 


XMMES1.374 


8.9731188 


-43.353933 


22.23 


20.56 


19.47 


16.07 


17.41 


0.25-0.80 


XMMES 1_444 


9.1195568 


-43.475508 


21.07 


19.44 


18.24 


14.63 


16.30 


0.30-0.80 



1 Magnitude in the Johson-Cousins Vega system. 





Fig. 6. Radial distribution of the background subtracted surface brightness, normalized to the peak central value, 
of the sources XMMES1_145, XMMES1.148, XMMES1.184, XMMES1.224, XMMES 1.363, XMMES1.374 and 
XMMES1.444, computed in concentric annuli (full dots), compared that of point-like sources at the same off-axis 
angles (open dots). 



the radial profiles of the point-like sources are significant 
at more than 99.9%. 

Figure \5\ shows the X-ray contours of the 7 sources 
overlayed on deep R band optical images obtained using 
VIMOS@VLT ESO telescope. We see that all sources can 
be identified with clusters or groups of galaxies. 



We extracted the X-ray spectra of the 7 extended 
sources from circular regions with radii in the range 35" 
to 70". The background counts were extracted from the 
nearby source- free regions. The response and ancillary files 
were generated by the XMM-SAS tasks, rmfgen and ar- 
fgen, respectively. The MOS1 and MOS2 spectra are 
combined together. The spectra were fitted using XSPEC 
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Fig. 7. PN and MOS counts spectra of the two brighest extended sources, XMMES1_145 (left panel) and XMMES1_224 
(right panel), along with their best fit thermal models. 




XMMES1_224 



XMMES' 



45 



2.5 



3.5 



kT 



keV 




Fig. 8. Left panel: the 68%, 90% and 99% redshift-temperature \ 2 contours for the two clusters of galaxies 
XMMES1_145 and XMMES1_224. Right panel: The 68%, 90% and 99% redshift-metal abundance X 2 contours. 



(version 11.3.1. lArnaud 2003(1 . to a thin plasma thermal 
model with free temperature and redshift. For the two 
brightest sources (XMMES1J.45 and XMMES1_224, see 
Table 4), the statistics is good enough to let also the metal 
abundance to be free to vary in the fit. Figure [7] shows the 
PN and MOS spectra of these two sources along with the 
best fit thermal model. In the other five cases, we fixed 
the metal abundance to the typical value (0.3) found in 
clusters of galaxies over broad redshift range (see, e.g., 
ITozzi et al. "20 031. Table 4 gives for each of the extended 
sources the X-ray position, the net counts in the 0.5-10 
keV energy range and the corresponding flux, the best fit 
temperature, redshift and X-ray luminosity. 



Figure [3] shows the temperature-redshift and redshift- 
metal abundance \ 2 confidence contours for the sources 
XMMES1_145 and XMMES1.224. The quality of the 
XMM-Newton data is such to tightly constrain the red- 
shift of the two bright clusters (0.21±0.01 and 0.39±0.02). 
Figure [5] shows the temperature-redshift x 2 contours for 
the sources XMMES1_374, XMMES1_363, XMMES1_184 
and XMMES1_148. Although in these cases the statis- 
tics is limited (see counts and fluxes in Table 4) the 
XMM-Newton data are still good enough to provide 
rough information on the redshift of the sources. Note 
that the redshift of XMMES1_363 is constrained to be 
> 0.7 at the 90% confidence level, and that 4 of the 
sources ( XMMES1.224, XMMES1.374, XMMES1.184 
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kT keV 



Fig. 9. The 68%, 90% and 99% redshift-temperature X 2 
contours for the sources XMMES1_148, XMMESl_184, 
XMMES1.363 and XMMES1.374. 



already noted by Rosati et al. (2002). The ELAIS-S1 de- 
termination is consistent within the errors with those of 
Moretti et al. (2003) and Baldi et al. (2002). The ELAIS- 
Sl 2-10 keV number counts are compared to those of 
Rosati et al. (2002), Moretti et al. (2003), Baldi et al. 
(2002) and Giommi et al. (2000) (BeppoSAX). All deter- 
minations are consistent with each others at fluxes lower 
than 2 - 3 x 10~ 14 erg cm~ 2 s~ x . The ELAIS-S1 5-10 
keV number counts are compared to the those of Rosati et 
al. (2002), Hasinger et al. (2001) (Lockman Hole), Baldi 
et al. (2002) and Fiore et al. (2001) (BeppoSAX). The 
ELAIS-S1 number counts are slightly lower than the other 
determinations but still consistent within the errors. In 
summary, the agreement of the ELAIS-S1 number counts 
with other determinations is reasonably good, implying 
that the ELAIS-S1 region, at least from this point of view, 
can be considered as representative of the average X-ray 
source population. 



and XMMES1_148) are consistent with having the same 
redshift (z~0.4). The luminosities and temperatures in 
Table 4, both determined through the X-ray spectral 
fittings, are in good agreement with the luminosity- 
temperature relationship of groups and clusters of galaxies 
(see, e.g.. lArnaud 20051 and references therein). 

Table 5 gives the position and the B, V, R, K and 
J magnitudes of the brightest galaxies of each of the 
seven X-ray sources in Figure |3J We estimated their pho- 
tometric redshifts by comparing the observed optical to 
near infrared spectral energy distribution with the Rocca- 
Volmerage galaxy templates (see, e.g., Fioc & Rocca- 
Volmerage 1997). Unfortunately, the available photometry 
does not allow a precise determination of the redshifts of 
the galaxies, but the redshift confidence intervals in Table 
5 are in all cases consistent with the redshift estimates in 
Table 4. In particular, the photometric redshifts confirm 
that XMMES1_363 is at a redshift higher than the rest of 
the extended sources sample. 

3. The X-ray number counts 

The XMM-Newton survey in the ELAIS-S1 area covers a 
range of fluxes over an area which is large enough to con- 
strain source counts in a flux range bridging deep, pencil 
beam, and shallower, large area surveys. Figure ITU1 plots 
the integral number counts in the F, S, H and HH bands. 
Table 6 gives the number counts in the S, H and HH 
bands, which can be more easily compared with previ- 
ous determinations, together with the sky coverage in the 
same bands. Previous determinations of the S, H and HH 
number counts are over-plotted on the ELAIS-S1 number 
counts. The 0.5-2 keV number counts are compared to 
those of Rosati et al. (2002) (Chandra Deep Field South, 
CDFS), Moretti et al. (2003) (compilation), Baldi et al. 
(2002) (HELLAS2XMM) and Hasinger et al. (2005) (com- 
pilation). The CDFS number counts are the lowest, as 



Table 6. Integral number counts and sky coverage 



Flux 

1(T 14 erg cm" 2 s _1 


N 


Counts 
deg" 2 


Sky-coverage 
deg 2 


0.5-2 keV 


0.10 


23 


854±60 


0.03 


0.15 


50 


564±31 


0.38 


0.23 


61 


473±27 


0.57 


0.36 


84 


321±22 


0.62 


0.55 


59 


223±19 


0.63 


0.83 


48 


138±15 


0.64 


1.30 


30 


75±11 


0.64 


1.99 


23 


33±7 


0.64 


3.04 


5 


11±4 


0.64 


4.55 


1 


5±3 


0.64 


7.15 


1 


3±2 


0.64 


11.00 


1 


2±2 


0.64 


2-10 keV 


0.47 


18 


746±91 


0.04 


0.63 


19 


470±43 


0.10 


0.86 


26 


297±25 


0.22 


1.16 


32 


212±19 


0.42 


1.57 


35 


152±16 


0.54 


2.12 


27 


99±13 


0.59 


2.87 


22 


46±9 


0.61 


3.88 


15 


21±6 


0.61 


5.25 


1 


5±3 


0.61 


7.10 


1 


2±2 


0.61 


5-10 keV 


0.69 


3 


156±47 


0.013 


0.96 


1 


91±24 


0.07 


1.34 


3 


51±12 


0.17 


1.86 


7 


35±9 


0.31 


2.59 


11 


15±5 


0.45 


3.60 


1 


2±2 


0.52 



Detected sources per flux bin. Note that a few sources fall 
out of the flux range, which was used to evalued the integral 
number counts in each energy band. 
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Fig. 10. The source number counts in the 0.5-10 keV F band (upper left panel), 0.5-2 keV S band (upper right panel), 
2-10 keV H band (lower left panel), 5-10 keV HH band (lower right panel). The Figure also shows a few number counts 
from literature: the 0.5-2 keV number counts from Rosati et al. (2002), Moretti et al. (2003), Baldi et al. (2002) and 
Hasinger et al. 2005 (upper right panel); the 2-10 keV number from Rosati et al. (2002), Moretti et al. (2003), Baldi 
et al. (2002) and Giommi et al. (2000) (lower left panel); the 5-10 keV number from Rosati et al. (2002), Hasinger et 
al. (2001), Baldi et al. (2002) and Fiore et al. (2001) (lower right panel). 



4. X-ray spectral properties 

The broad band X-ray spectral properties of the ELAIS-S1 
sources can be investigated using hardness ratios. Many of 
the ELAIS-S1 sources have a high hardness ratio HR=(H- 
S)/(H+S) (where H and S are the fluxes in the hard and 
soft band, respectively, (see, e.g., Fiore et al. 2003)), in- 
dicating a hard spectrum. H-S/H+S for the full 2-10 keV 
sample, 205 sources, is plotted in Figure ^] left panel, 
as a function of the 2-10 keV flux. Lines of constant Njj 
are overlaid. Following Hasinger et al. (2003) a threshold 
of the (H-S)/(H+S) ratio can be used to roughly separate 
obscured from unobscured sources. To this purpose we use 
HR=0.5, which corresponds to log N# > 21.5, 22.2, 22.7 
at z=0, 1, 2, respectively (assuming a power law slope 
a=0.8). 

According to this definition, the fraction of obscured 
sources to the total is 0.36. Of course, since errors on HR 
are large, especially at low fluxes, the separation between 
obscured and unobscured sources has only a statistical 
meaning. To evaluate how the fraction of the obscured 



sources is affected by the errors of HR, we generated 
1000 random source distributions with the same number 
of sources as the real data, by assigning to each source 
a random HR value drawn from a gaussian distribution 
with mean equal to the measured HR and a equal the HR 
error. We then computed the dispersion of the number of 
sources with HR> 0.5, which turned out to be about 5%. 
The dotted histogram in the right panel of Figure ^] is 
the average of the simulated (H-S)/(H+S) distributions. 
The last point of the simulated (H-S)/(H+S) distribution 
deviates from that of the real (H-S)/(H+S) distribution 
because the former takes properly into account the upper 
limits, which are plotted in the solid histogram at their 
nominal values. 

Adopting the above definition, Figure Il2l shows num- 
ber counts of obscured and unobscured sources respec- 
tively, while Figure El skews the ratios between the num- 
ber of the obscured sources to the total. The fraction of 
obscured sources increases with the energy of the X-ray 
band, consistently with Comastri et al. (2001). In the 2- 
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Fig. 11. Left panel: The hardness ratio (H-S)/(H+S) as a function of the flux for the whole 2-10 keV sample. Typical 
hardness ratio errors are shown at different fluxes, la upper limits on the S band flux are shown with by arrows. 
Large filled circles represent the median and interquartile in two flux ranges, evalued using the so called "survival 
analysis" ( Ka plan and Meier, 1958[ |Miller 1981, p. 74| ) . They are compared with the same determination from the 
HELLAS2XMM sample (open squares). Constant lines are the hardness ratios expected from a power law model with 
slope a =0.8, reduced at low energy by photoelectric absorption from neutral gas (Nff= 10 22 cm -2 and 10 23 cm -2 ) at 
three representative redshifts: z=0, z=l and z=2. Right panel: The distribution of (H-S)/(H+S) for the same sample 
in the left panel (solid histogram). Upper limits are identified by the shaded histogram. The dotted histogram is the 
average of the simulated (H-S)/(H+S) distributions. 



10 keV band the fraction of obscured sources increases 
from ^20% at fluxes > (4 - 5) x 10~ 14 erg cm~ 2 s _1 
to ~ 40% at fl uxes ~ (1 - 2) x 10~ 14 erg cm" 2 s" 1 
(consistent with IPiconcelli et al. 20031 lUeda et al. 2 003 
IPerola et al. 20041 ILa Franca et al. 2005|l . It should be 
beared in mind that all this refers to moderately obscured 
sources only, since sources obscured by column densities 
N# >a few 10 24 cm -2 , the so called Compton thick AGNs, 
are rarely detected even in 2-10 keV samples (see Tozzi et 
al. 2006). 



5. The angular correlation function 

We now focus on the clustering properties of our X-ray 
sources. The fraction of active galaxies in our samples is 
probably higher than 80-90% in all bands (according to 
previous identification of X-ray sources at similar flux lim- 
its, see, e.g., Mainieri et al. 2002, Fiore et al. 2003, Barger 
et al. 2005). We computed the two point source angular 
correlation function fACF. IPeebles 1980) in the three S, H 
and F bands using the Landy & Szalay (LS, 1993) estima- 
tor. The ACF function is given by the following relation: 

... fci x DD -2 x k 2 x DR + RR 



where DD is the number of pairs of sources in the 
real data with separation angle in the interval 9,9 + d9, 
RR is the number of pairs in the same separation angle 
interval using a spatially random distribution of sources, 
DR is the number of pairs obtained comparing real data 
and random data, and fci = " r (" r — ii fc 2 = i) w j iere 

' 1 n d (n d -l)' z 2xn d ' 

n^ is the total number of real sources and n r is the number 
of random sources. 

To evaluate w(9) we followed two different methods: 
(1) we generated 300 random source distributions with 
the same number of sources as the real data by convolving 
a random source distribution with the differential source 
number counts and with the appropriate sensitivity maps. 
We adopted a cut in flux for both the real source and ran- 
dom source distributions to exclude from the ACF anal- 
ysis the 10% faintest sources in the samples. We then 
computed the ACFs according to eq. (1) for each ran- 
dom realization and in logarithmically spaced 9 bins. The 
distribution of the 300 ACF values in each 8 bin is in 
most cases reasonably well fitted by a gaussian distribu- 
tion and therefore we computed the mean and the rms 
values in each 9 bin. (2) we generated a single random 
source distribution using 100 times more sources than in 
the real samples. Also in this case, we adopted a cut in 
flux for both the real source and random source distribu- 
tions to exclude from the ACF analysis the 10% faintest 
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Fig. 12. The source number counts in the F band (upper 
left panel), S band (upper right panel), H band (lower left 
panel) , HH band (lower right panel) . The open circles indi- 
cate the source number counts of the unobscured sources. 
The solid circles indicate the source number counts of the 
obscured sources, as defined in Sect. 4. 
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Fig. 13. Ratio between the number of the obscured 
sources to the number of the total sources in the S, H, 
HH bands (from top to bottom). 



sources and computed the ACFs according to eq. (1). In 
this case errors were computed adopting the following for- 
mula (from Landy & Szalay, 1993): 



6(6) 



(l + w(9)) 



(2) 



ViDD) 

Figure PHI shows the ACF in the S and H bands com- 
puted in bins of log 9 = 0.075 following the second 
method. The solid lines represent the best fit w(6) = 
(d/Oo) 1 1 — ujq. function to the ACF. The ujq factor has 
been introduced following Landy & Szalay (1993), Daddi 
et al. (2000) and Basikalos et al. (2004) to take into ac- 
count the so called "integral constraint" . This results from 
the fact that the correlation function is evaluated in a lim- 
ited area, while the ACF is normalized to zero over the full 
sky. According to Roche et al. (1999)]this implies that the 
w(9) computed using eq. (1) over-estimates the real w(6) 
by a constant quantity uq= (l/flp) 1 " 7 x C, where C has 



been estimated following ag ain Roc he et al. (1999) Fixing 
1.8, as in Vikhlinin & Forman (1995) Basilakos et 



7 to 

al. (2005) and D'Elia et al. (2005), produces the follow- 
ing best fit correlation angles: 6q = 5.2 ± 3.8 arcsec, 
6» = 12.8 ± 7.8 and 6 = 7.5 ± 3.1 arcsec in the S, H and 
F bands respectively (errors are at la confidence level). 
Letting 7 to vary of course increases the already large er- 
rors on Qq but does not change their best fit values. In fact 
the best fit 7 turn out to be and 1.8±g*a for the 

S and H bands respectively. To test the robustness of our 
determinations we computed the ACF with different bin 
sizes from log 9 = 0.05 to log 9 = 0.1. For the 0.5-2 keV 
ACF we found best fit 9q in the range ~4-9 arcsec, while 
for the 2-10 keV ACF we found best fit values in the range 
~9-19 arcsec. In all cases the best fit #0 were consistent, 
within the errors, with the values given above. 

Very similar results are obtained by using the first 
method to compute the w(9) and the corresponding er- 
rors. 

Since the correlation angles we probe are of the same 
order of magnitude of the XMM-Newton PSF, the so 
called amplification bias (see, e.g., Vikhlinin & Forman 
1995) could in principle affect our estimates of w(9) at 
small separations, due to sources closer than a few arc- 
sec being observed as a single source. However, we note 
that at our limiting flux the source confusion affects our 
images little (see Sect. 2.5). Furthermore, since our fits 
have been performed on angular scales 8 > 100", much 
greater than the XMM-Newton PSF, the amplification 
bias is not expected to affect our estimates of the ACF. 
Indeed, Basilakos et al. (2005) found that amplification 
bias is negligible in their XMM-Newton 2dF survey. 

5.1. Source clustering 

The ACF computed in the previous Sect, measures the 
strength of the angular correlation, but does not say any- 
thing about where the sources are clustered. To investi- 
gate this issue we plot in Figure H3k') and ED 3 ) the source 
density map of the 0.5-2 keV and 2-10 keV sources, respec- 
tively, smoothed with a box of 10 arcmin side. The high 
density regions of both maps encompass the two clusters 
at z~ 0.4 (XMMES1.224 and XMMES1.374) and other 



Puccetti et al.: The XMM-Newton survey of the ELAIS-S1 field 



15 



respectively using the following equation: 
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Fig. 14. The angular correlation function in the S band 
(top panel), and H band (bottom panel) computed using 
the Landy & Szalay estimator. The solid lines represent 
our best fit to the data between 100 and 1000 arcsec. 



two clusters or groups of galaxies with a redshift consistent 
with 0.4. However, these observed source density maps are 
affected by instrumental biases, i.e. a higher source den- 
sity is visible in the regions with the highest exposure 
(see Figure EJl . To account for this effect, we computed 
the deviations of the real source distribution from a ran- 
dom source one, which has been obtained by convolving 
a random source distribution with the differential source 
number counts and with the appropriate sensitivity maps. 
We plot in Figure IT31 c) and 1151 d) the normalized source 
density NNSD of the 0.5-2 keV and 2-10 keV sources 



NNSD(X, Y) = (NSD(X, Y) ~ NRD(X, Y))/NRD(X, Y) 

(11) 

where NSD is the real source density and NRD is the 
density of a random source sample. The resulting normal- 
ized source density map has then been smoothed with a 
box of 10 arcmin side. The high density regions of the 
0.5-2 keV map still encompass the four extended sources 
clusters at z~ 0.4. On the other hand, the highest density 
region of the 2-10 keV map is just located in the central- 
west region of the field, in an area which coincides with 
one of the high density regions of the 0.5-2 keV map, but 
in which no extended source has been found. 



6. Discussion 

We have surveyed with XMM-Newton the central ~ 0.6 
deg 2 region of the ELAIS-S1 field down to flux limits of 
~ 5.5 x 10~ 16 erg cm" 2 s" 1 (0.5-2 keV), - 2 x 10~ 15 
erg cm~ 2 s _1 (2-10 keV) and ~4x 10~ 15 erg cm~ 2 s _1 
(5-10 keV ). 448 sources have been detected in the 0.5-10 
keV F band images with a probability lower than 2xl0 -5 
that they are Poisson fluctuations of the background. 395, 
205 and 31 sources are detected in the 0.5-2 keV, 2-10 keV 
and 5-10 keV band images respectively. Of these, 26 and 4 
sources are detected only in the S and H band but not in 
the F band. The total number of detected sources is 478. 

The source number counts are in good agreement with 
previous determinations from large area surveys and com- 
pilations (Baldi ct al. 2002, Moretti et al. 2003), show- 
ing that field to field variance is smoothed out at the 
scale of the ELAIS-S1 survey (0.6 deg 2 ) and implying 
that the ELAIS-S1 region, at least from this point of 
view, can be considered as representative of the aver- 
age X-ray source population. The average spectral prop- 
erties of the ELAIS-S1 sources are consistent with those 
of other surveys covering similar or larger areas at a sim- 
ilar depth. In particular, we find that: a) the fraction of 
obscured sources increases with the energy of the X-ray 
band and b) in the 2-10 keV band this fraction increases 
from ~20% at fluxes > (4 - 5) x 10~ 14 erg cm -2 s _1 
to - 40% at flu xes ~ (1 - 2) x IP" 14 erg cm" 2 s" 1 
(consistent with IPiconcelli et al. 20031 lUeda et al. 20031 
IPerola et al. 20041 ITa Franca et al. 2005jl . 



We identified 7 clearly extended sources in the field and 
estimated their redshift through spectral fits with ther- 
mal models. In two cases the constraint on the redshift is 
quite accurate, <j(z)/(1 + z) < 1.5%, in three other cases 
a(z)/(l + z) ~ (4 — 10)%, and in the remaining two cases 
the upper limit to the redshift is not well determined. 
Of course these uncertainties do not take into account 
possible systematic errors and should therefore be con- 
sidered with some caution. The redshift of XMMES1_363 
is constrained to be > 0.7 at the 90% confidence level. In 
four cases (XMMES1.224, XMMES1.374, XMMES1.184 
and XMMES1_148) the redshift is consistent with z~0.4. 
Photometric redshifts estimated using the B, V, R, K and 
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Fig. 15. Top left panel a): the 0.5-2 keV source density contours, corresponding to -30%, -15%, 5%, 20% and 40% 
of the average source density; top right panel b): the 2-10 keV source density contours corresponding to -40%, -20%, 
5%, 30% and 55% of the average source density. Darker regions correspond to undcrdensities and lighter regions to 
overdensities. The diamonds in all panels represent the 7 extended sources in the field (see Sect. 2.5). Bottom left 
panel c): same as top left panel for the normalized 0.5-2 keV source density; bottom right panel d): same as top right 
panel for the 2-10 keV source density contours. 



J magnitudes of the brightest galaxies associated with the 
seven X-ray sources are all consistent with the redshifts 
estimated through their X-ray spectra. Interestingly, a 
peak at z^ 0.4 is present in the spectroscopic redshift 
distribution of radio sources around XMMES1_374 and 
XMMES1_184 (Gruppioni ct al. in preparation). The 0.5- 
10 keV luminosity of XMMES1.224 is - 8 x 10 43 erg s" 1 , 
while the luminosity of the other three extended sources 
at a similar redshift is 15-60 times smaller. The angular 
separations between the luminous cluster XMMES1.224 



and the 3 lower luminosity clusters or groups of galax- 
ies are 5' (XMMES1.184), 7.5' (XMMESU48) and 15.5' 
(XMMES1.374), corresponding to 1.6, 2.4 and 5 Mpc at 
z~0.4, for the assumed cosmology. A detailed analysis of 
the extended sources and of the large scale structure at 
z~ 0.4, including spectroscopic redshifts, will be presented 
in a forthcoming publication. 

We have computed the angular correlation function of 
the sources in the S and H bands using the Landy & Szalay 
(1993) estimator. Fixing the slope of the correlation func- 
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tion at 7 =1.8, we found the following correlation angles: 
#o = 5.2 ± 3.8 arcsec and #o = 12.8 ± 7.8 arcsec in the S, 
and H bands respectively, while the correlation angle in 
the 0.5-10 keV band is intermediate between these values. 
These estimates are fully consistent with the statistically 
poorer D'Elia et al. (2005) determinations. Basilakos et al. 
(2005) and Basilakos et al. (2004) found 6 = 10.4 ± 1.9 
arcsec in the S band and 6o = 28 ± 9 in the H band for 
samples of sources of size similar to our ELAIS-S1 sample, 
but with a flux limit — 5 — 6 times higher than ours and 
spread over an area 4 times wider. 

The correlation angle of hard X-ray selected sources 
is formally —2.5 times larger than that of soft X-ray se- 
lected sources (a correlation angle of hard sources higher 
than soft sources has been suggested also by Yang et al. 
2003 and Basilakos et al 2004), although the difference is 
significant at only ~ la level. If real, this difference may be 
due to the different redshift distributions of hard and soft 
sources, because the angular correlation angle depends 
also on the average distance between galaxies, and there- 
fore on their redshift. For a proper comparison one would 
need to evaluate the spatial correlation function, once the 
redshifts of the sources are known. A campaign to ob- 
tain optical spectroscopy of the optical counterparts of the 
X-ray sources is currently on going using VIMOS@VLT 
and the results will reported in a future publication. For 
the time being, we can obtain a rough estimate of the 
present-day correlation length ro assuming an appropri- 
ate redshift distribution dN/dz and inverting the Limber 
equation IjLimber 19531 [Pee bles 1980). We have estimated 
the observed redshift distributions of the S and H band 
sources, by convolving the best fit luminosity functions of 
La Franca et al. (2005, 2-10 keV) and Hasinger et al. (2005, 
0.5-2 keV) with ELAIS-S1 sky coverage. Of course this ap- 
proach does not take into account the presence of signif- 
icant structures at a given redshift in the ELAIS-S1 field 
which could bias the redshift distribution of the ELAIS-S1 
sources (one such structure may be present at z— 0.4, see 
Sect 2.5), and therefore the following results are should 
be considered as upper limits. Assuming an evolutionary 
parameter e 1 = 7 — 3 = —1.2 (for 7 = 1.8, the so called 
"comoving clustering model" ) , which should be appropri- 
ate for active galaxies IfKundic 19971 ICroom et al. 20(711 
lUroom et al. 2005jl . we obtain r = 12. 8tf® h" 1 Mpc and 
ro = 17.9^g'i n_1 Mpc in the S and H band respectively. 
So far we have obtained secure optical spectroscopy red- 
shifts for 87 and 49 sources in the S and H bands (Feruglio 
et al., in preparation). Using these (admitedly incomplete) 
redshift distributions, which clearly shows strong peaks at 
z=0.4, we obtain r = 9. 8t.ll h_1 M P C and r o = 13 - 4 -l6 
h _1 Mpc in the S and H band, respectively. 

The S and H band redshift distributions are peaked at 
z~l and z— 0.9 respectively (median redshift) and there- 
fore we are effectively measuring the clustering at those 
redshifts. The correlation length at a given redshift is: 



and therefore ro(z) = ro for 



ro(z) = r (1 + z) i 
e = -1.2 and 7 = 1.8. 

Our correlation lengths are lower but still consistent 
within the large errors with those of Basilakos et al. (2005) 
and Basilakos et al. (2004). On the other hand, for 7 = 1.8 
Gilli et al (2005)] find r -12 fr 1 Mpc and r 



^6 1T 



1 e parameterizes the redshift evolution of the clustering (see, 
e.g., Overzier et al. 2003). 



Mpc for the faint sources in the CDFS and CDFN fields 
respectively. The former estimate is consistent with our 
determination. We also note that although the median 
redshift of these sources is similar to the median redshift 
assumed above for the ELAIS-S1 sources, their fluxes (and 
therefore luminosities) are a factor of 5-10 lower than 
those of the sources in the ELAIS-S1 field. Yang et al. 
(2006) present a spatial correlation function analysis of 
the non-stellar X-ray point sources in the Chandra Large 
Area Synoptic X-ray Survey of Lockman Hole Northwest 
(CLASXS), covering a total area of 0.4 deg 2 down to a 
flux limit of 3xl0~ 15 erg cm~ 2 s _1 in the 2-8 keV band. 
They find r = 8.IIV2 n_1 Mpc, consistent with our de- 
termination obtained using the observed, but incomplete, 
redshift distributions. 

We have also compared our correlation lengths to those 
obtained from optically selected QSOs, optically and near 
infrared selected galaxies and radio sources. Croom et 
al. (2005) find a low correlation length in redshift space 
so = 4.8 ± 1 fr 1 for 7 = 1.23. Grazian et al. (2004) find a 
low ro = 8.6 ± 2.0 h _1 Mpc for 7 = 1.56 and for a sample 
of local (0.02<z<0.22), optically selected AGN. Both de- 
termination are lower than our values. On the other hand, 
the correlation lengths of the ELAIS-S1 X-ray sources are 
comparable with that of Extremely Red Objects (EROs) 
(see, e.g., Daddi et al. 2001, r = 12 ± 3 hr 1 Mpc for 
7 = 1.8 with l<z<1.2) and luminous radio sources (see, 
e.g., Overzier et al. 2003, r = 14 ±3 h" 1 Mpc for 7 = 1.8 
) at similar redshifts (z— 1). 

We finally investigated whether the source density and 
clustering (on scales between 3 and 6 arcmin, where the 
signal in the ACF is more significant), are spatially related 
to the large scale structure discovered at z— 0.4. Figure 
IT5l shows that there is a high normalized source density 
for the S band sample around the clusters and groups at 
z— 0.4 (or with a redshift consistent with this value) ex- 
tending both toward East and toward South/ West, sug- 
gesting that the structure is complex, with a size com- 
parable to the full XMM-Newton coverage. The highest 
normalized source densities of the H band sources are lo- 
cated in the South West and central- West regions of the 
field. The first is nearly coincident with one of the peaks 
of the S band source density, while the second has not a 
S band counterpart. The H band density peaks appear to 
be unrelated to the position of the extended sources at 
z~ 0.4. 

Sensitive observations over a wider area are clearly 
needed to both improve the statistics of the correlation 
functions and to cover several structures of the size of 10- 
20 Mpc in the same area, thus allowing us to understand 
whether these are peculiar structures, found by chance in 
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selected areas of the sky, or whether, rather, they are rep- 
resentative of the average Large Scale Structure. The on- 
going COSMOS multiwavelength survey will achieve these 
goals in the following years. 

Acknowledgements. Part of this work was supported by 
MIUR COFIN-03-02-23 and INAF/PRIN 270/2003. We thank 
Adriano Fontana for allowing us to use his photometric redshift 
code. 

References 

Arnaud M. Proc. Enrico Fermi, International School of 

Physics Course CLIX, eds. F. Melchiorri & Y. Rephaeli, 

astro-ph/0508159 
Ar naud K. A. 2003, 'XSPEC User Guide for version 11.3.x ' 

http: / / xspec.gsfc.nasa.gov / docs /xanadu /xspec /xspecl 1 / 

manual / manual.html 
Baldi, A., Molendi, S., Comastri A., et al. 2002, ApJ, 564, 190 
Barger A., Cowie L., Mushotzky, R. F. et al. 2005, AJ, 129, 

578 

Basilakos, S., Plionis, M., Georgakakis, A., Georgantopoulos, 

I. 2005, MNRAS, 356, 183 
Basilakos, S., Georgakakis, A., Plionis, M., Georgantopoulos, 

I. 2004, ApJ, 607, 79 
Berta et al. 2006, submitted to a 
Cavaliere, A. & Vittorini, V 2000, ApJ, 543, 599 
Comastri, A., Fiore, F., Vignali, C. et al. 2001, MNRAS, 327, 

781 

Croom, S. M.; Boyle, B. J.; Shanks, T. et al. 2005, MNRAS, 
356, 415 

Croom, S. M., Shanks, T., Boyle, B. J. et al. 2001, MNRAS, 
325, 483 

Daddi, E., Broadhurst, T., Zamorani, G. et al. 2001,A&A, 376, 
825 

Daddi, E., Chnatti, A., Pozzetti, L. et al. 2000, A&A, 361, 535 
Damiani, F.; Maggio, A.; Micela, G., Sciortino, S., 1997, ApJ, 
483, 350 

D'Elia, V.; Fiore, F.; Elvis M. et al. 2004, A&A422, 11 
Fioc, M. & Rocca-Volmerange, B. 1997 A&A 326, 950 
Fiore, F., Brusa, M., Cocchia, F. et al. 2003, A&A, 409, 79 
Fiore, F., Giommi P., Vignali C. et al. 2001, MNRAS, 327, 771 
Franceschini, A., Hasinger, G., Miyaji, T., &Malguori, D. 1999, 

MNRAS, 310, L5 
Gilli, R., Daddi, E., Zamorani, G. et al. 2005, A&A430, 811 
Giommi, P., Perri, M., Fiore, F. 2000, A&A, 362, 799 
Ghizzardi, S. "In-flight calibration of the PSF for the PN cam- 
era", EPIC-MCT-TN-012, 2002 
Ghizzardi, S. "In-flight calibration of the on-axis and near off- 
axis PSF for the Mos-1 and Mos-2 cameras", EPIC-MCT- 
TN-011, 2001 

Grazian, A., Negrello, M., Moscardini, L. et al. 2004 AJ, 127, 
592 

Groth, E. J. & Peebles, P. J. E. 1977, ApJ, 217, 385 
Gruppioni, C, Ciliegi, P., Rowan-Robinson, M. et al. 1999, 

MNRAS, 305, 297 
Hasinger, G., Miyaji, T., Schmidt, M. 2005 A&A, 441, 417 
Hasinger, G. 2003, AIP Conf. Proc. 666,227 
Hasinger, G.; Altieri, B.; Arnaud, M. et al. 2001, A&A, 365, 

L45 

Huchra J. P., Geller M. J., 1982, ApJ, 257, 423 
Kaplan E. L. and Meier P., 1958, J. A. Statistical Association, 
53, 457 



Kundic, T. 1997, ApJ, 482, 631 

La Franca F., Fiore F., Comastri A. et al. 2005, Scheduled on 
ApJ, 635, 864 

La Franca F, Gruppioni, C, Matute, I. et al. 2004, AJ, 127, 
3075 

Landy, S.D., Szalay, S.D. 1993, ApJ, 412, 64 
Lari, C, Pozzi, F., Gruppioni, C. et al. 2001, MNRAS, 325, 
1173 

Limber, D. N. 1953, ApJ, 117, 134 

Lonsdale, C; Polletta, M. d. C; Surace, J. et al. 2004, ApJS, 
154, 54 

Lonsdale, C; Smith, H. E.; Rowan- Robinson, M. et al. 2003, 
PASP, 115, 897 

Lumb, D. H.; Finoguenov, A.; Saxton, R. et al. 2003, ExA, 15, 
89 

Mainieri, V., Bergeron, J., Hasinger, G., Lehmann, I., Rosati, 
P., Schmidt, M., Szokoly, G., & Delia Ceca, R. 2002, A&A, 
393, 425 

Menci N., Fiore, F., Perola, G.C, Cavaliere, A. 2004, ApJ, 606. 
58 

Miller R.G.Jr. 1981, Survival Analysis (New York:Wiley) 
Moretti A., Campana S., Lazzati D., Tagliaferri, G. 2003, ApJ, 
588, 703 

Oliver, S.; Serjeant, S.; Efstathiou, A. et al. 2000, LNP, 548, 
280 

Overzier, R. A., Rttgering, H. J. A., Rengelink, R. B., Wilman, 

R. J. 2003, A&A, 405, 530 
Peebles P.J.E., 1980, The Large-Scale Structure of the Universe 

(Princeton: Princeton University Press) 
Perola, G. C, Puccetti, S., Fiore, F. et al. 2004, A&A, 421, 491 
Piconcelli, E., Cappi, M., Bassani, L.,Di Cocco, G., Dadina, 

M. 2003, A&A, 412, 689 
Pillitteri I., Micela G., Damiani F., Sciortino S. 2006, A&A, 

450, 993 

Read, A. M.; Ponman, T. J. 2003, A&A, 409, 395 

Roche N., Eales S. A., Hippelein H. and Willott C. J. 1999, 

MNRAS, 306, 538 
Rosati, P.; Tozzi, P.; Giacconi, R. et al. 2002, ApJ, 566, 667 
Rowan-Robinson, M.; Lari, O; Perez-Fournon, I. et al. 2004, 

MNRAS, 351, 1290 
Springel, V., White, S. D. M., Jenkins, A. et al. 2005, Nature, 

435, 629 

Struder, L.; Briel, U.; Dennerl, K. et al. 2001, A&A, 365, L18 
Tozzi, P., Rosati, P., Ettori, S. 2003 ApJ, 593, 705 
Tozzi, P., Gilli, R., Mainieri, V. 2006 A&A, 451, 457 
Turner, M. J. L.; Abbey, A.; Arnaud, M. et al. 2001,A&A, 365, 
L27 

Ueda Y., Akiyama M., Ohta K., Miyaji T. 2003, ApJ, 598, 886 
Vikhlinin A., & Forman W. 1995, ApJ, 455, L109 
Wolf, O, Meisenheimer, K., Kleinheinrich, M. et al. 2004, 
A&A, 421, 913 

Yang, Y., Mushotzky, R. F., Barger, A. J., L. L. Cowie 2006, 

submitted to ApJ, astro-ph/0601634 
Yang, Y., Mushotzky, R. F., Barger, A. J., et al. 2003, A&A, 

585, L85 



